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ABSTRACT 



Observations of the temporal evolution of the wave-formed ripples are analyzed 
in terms of geometric properties, migration rate, and forcing by wave and current 
velocities. Three weeks of bedform observations were obtained using imderwater video 
of a sheet of laser light projected on the bed at the Monterey Bay Inner Shelf Observatory 
in thirteen-meter depth water. The bed consists of fine sand (mean grain size 0.12mm). 
Low to moderate narrow-band swell waves occurred during the observation period. 
Ripple geometry consisted of orbital and suborbital vortex wave ripples, and relic ripples 
left after larger wave events. The bedforms generally changed size and shape when the 
grain roughness Shields parameter exceeded a critical value of 0.04 with the exception of 
the last event bringing into question the appropriate value of Ocru^n nature. Ripple 
migration was offshore at rates of 2-10 cm/day during active sediment transport events. 
Skevraess of velocities (low and high passed) were calculated to explain offshore ripple 
migration, but showed no preferred direction. Mean currents were weak and also showed 
no preferred direction. Significant positive (offshore) correlation was obtained between 
the short-wave envelope and inffagravity waves indicating wave-group forced bound long 
waves (surf beat) combined with stirring by the short waves might explain the offshore 
sediment transport and ripple migration. 
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I. INTRODUCTION 



As shoaling waves move from deep to shallow water over a sandy bottom, energy 
dissipation increases in the bottom boundary layer because of the presence of bed 
roughness in the form of wave-induced ripples. Ripple bedforms in coastal shallow water 
play a significant role in the interaction of bottom boimdary hydrodynamics and sediment 
transport. When wave energy is sufficient, fnctional stresses at the bottom cause wave- 
induced ripples to form. The flow over the ripples is dominated by the vortices created at 
the crest on the lee side of the ripples, which also act to enhance wave dissipation. 
Sediment transport is then dominated by avalanches of sand rolling down on the sides of 
the ripples, which is much larger than for a flat bed (Andersen 1999). Wave ripples, 
termed vortex ripples by Bagnold (1946), are triangular in shape with the slope of the 
sides being close to the angle of repose, usually shorter than 60 cm and not higher than 
about 6 cm (Engelund and Fredsoe, 1982). Ripples have complex geometry with 
different equilibrium shapes and sizes. 

Until recently, there have been few observations of bedform generated by waves 
and currents. Some of the first measurements of bedform generation due to waves and 
currents over a sandy bottom were by Inman and Bowen (1 963). Tanaka and Shuto 
(1984) performed laboratory experiments on bottom bedforms due to waves and currents 
interacting both in the cross-shore and along-shore directions, and were the first to 
provide an empirical formula for the prediction of ripple geometry due to waves and 
current forcing. 
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Most ripple wavelengths fall into three categories of wavelengths (Clifton and 
Dingier, 1 984). Orbital scale ripples, or “vortex” ripples, have wavelengths (X) that 
scale with wave orbital diameter (2A). Anorbital ripples have wavelengths that depend 
on sediment grain size (d) alone. Suborbital ripples represent a transitionary stage 
between orbital and anorbital ripples and depend on both wave orbital diameter and grain 
size. Important geometric properties to characterize ripples are their height (r|), steepness 
(ratio of ripple height to wavelength, 'r\IX), and cross and along shore symmetry. 
Symmetry in both these directions is an important property affecting formation and 
migration rates. A common problem with laboratory data is that the wave periods are 
generally short due to experimental apparatus limitations. Field studies have the 
limitations that the relationship between wave forcing and ripple geometry can be 
complicated due to the range of wave heights, periods, directions, and presence of low 
frequency ciirrents. 

Field experiments have been conducted using a variety of techniques to measure 
bedforms. Manual observations have been made by divers, which greatly limits the 
sample area and sample rate. Boyd et al. ( 1 988) and Amos et al. (1 988) used 
photographic measurements to investigate the role of waves and c;irrents in generating 
bedforms. Rotary sidescan sonar was used by Hay and Wilson (1994) and Traykovski et 
al. (1999) to measure the horizontal pattern of bedform evolution. Crawford and Hay 
(1998) developed an underwater video system consisting of a diode laser and a video 
camera to obtain two-dimensional structures of suspended sediment and bedforms. We 
have adapted a similar approach to analyzing ripple geometry and bedform evolution. 
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Measurements described here include waves, ripple formation, and increased 
wave attenuation. This study is a component of the ONR-sponsored shoaling waves 
experiment (SHOWEX) with the overall objective of measuring wave transformation 
mechanisms across the shelf An imderwater video camera in conjunction with a 
structured light source was used to obtain image sequences of the bottom, providing high- 
resolution digital pictures of ripples on the sandy bottom. The ability to resolve 1 mm- 
scale features in the bedforms gives this video imaging system the advantage over 
previous field experiments. The camera and light source were deployed in 13 meter depth 
at the Monterey Inner Shelf Observatory (Stanton 1999). 

Hypotheses for bed migration mechanisms on the inner-shelf (outside the 
surfzone) include wave skewness, mean cxirrents combined with wave stirring, and bound 
long waves. The non-linearity of shoaling waves causes waves to have an onshore 
skewness owing to more peaked crests and flatter troughs. This results in stronger 
onshore velocity and weaker offshore velocity with a net onshore sediment transport. 

Mean currents act to advect sediments that are stirred up and suspended by wave 
action. The direction of transport is in the direction of the current, which could be 
onshore, offshore, or alongshore. 

Boimd long waves act to transport sediments offshore. Natural waves are not 
regular, but have a variation in the wave height. A narrow-banded wave train consists of 
wave groups and associated bound long wave (Longuet-Higgins and Stewart 1964). Shi 
and Larsen (1984) recognized that bound long waves under wave groups could give an 
opposing contribution to the transport, because the bound long wave troughs occur under 
high waves Avith maximum offshore velocity coinciding "with high sediment 
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concentrations stirred by the wave group, while the onshore motion occurs under the low 
waves. Deigaard, Jakobsen, and Fredsoe (1999) found the nonlinearity of the turbulent 
boundary layer and the coupling between the offshore motion and high waves give an 
offshore mean bed shear stress. From their studies they concluded that bound long waves 
give an offshore contribution to sediment transport. 

The objective of this paper is to analyze observations of bedforms in relation to 
wave and current forcing. The paper is organized into five sections. The second section 
describes the experiment and the equipment that was deployed during September 2 to 
September 24, 1999, calibration and processing data that was acquired. The third section 
describes the results of ripple geometry and ripple migration analyses. A discussion of the 
results is included in the fourth section. The conclusion and recommendations finalize 
this thesis. 
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II. EXPERIMENT 



The MI SO instrument frame consists of a five-meter horizontal beam supported 
one meter above the bed by two vertical pipes jetted into the bed as displayed in Figure 1 . 
Mounted on the frame is a Broadband Acoustic Doppler Current Profiler (BADCP), 
under-water video camera, paro-scientific pressure sensor, scanning acoustic altimeter. 
Bistatic Coherent Doppler Velocity and Sediment Profiler (BCDVSP) and structured 
laser light. MISO is a component of the Rapid Environmental Assessment Laboratory 
(REAL) and was established July 26, 1999 and continues to be operational. MISO is 
located 400 meters off Del Monte Beach adjacent to the Naval Postgraduate School 
Oceanography laboratory. A local map depicting the location of MISO is shown in 
Figure 2a, and a regional map of California’s central coast is shown in Figure 2b. Waves 
incident to the MISO location in the Southern Monterey Bay tend to be narrow-banded in 
frequency and direction owing to the strong refraction as they pass over the Monterey 
Bay Submarine canyon, and protection by Point Pinos headland. The MISO array is 
located on a sand lens consisting of fine sand with shale outcrops 50 meters on either 
side. The sand at the MISO site was well-sorted quartz with a 0.125 mm median grain 
diameter (<^ 50 ) (Figure 3). 

The data analyzed here is for a 22-day period from September 2, to September 24, 
1999. This period provided comprehensive data from the underwater video camera, 
BADCP, and paro-scientific pressure sensor. The mean water depth was 12.8 m varying 
from 11.7mtol3.8m due to diurnal and semi-diurnal tidal forcing. No significant 
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storms occurred during this time period. The significant wave height peaked at 1 . 1 m and 
averaged 0.5 m. 

A. BROAD BAND ACOUSTIC DOPPLER CURRENT PROFILER (BADCP) 
Three-component velocity profiles with 1-m range resolution were measured 
using the BADCP. The instrument was situated approximately 1 meter off the seafloor 
with the center of the first range bin located 2.1 meters from the transducer face. The 
data were sampled at 0.976 Hz. Spectra were computed using thirty-minute ensembles. 
The spectra were averaged into fifty logarithmic-spaced frequency bins from 0.01 Hz to 
0.35 Hz. Mean-square wave bed velocities were calculated using a linear theory transfer 
function (/) to transform measured velocity at 4. 1 m, (/) , to the bed over a 0.05 
to 0.3 Hz frequency band 

( 1 ) 

.05 

where Hf , (/) = + z))’ ^ wavenumber and h+z is distance above the 

bed. The rms orbital excursion at the bed was computed from the band passed 
velocity spectrum at the bed using the integrating operator, (/) = ^ , on the 
frequency spectrum 

( 2 ) 

V0.05 

To determine the importance of mean currents on ripple formation, thirty-second 
cross-shore and alongshore mean currents were calculated from the BADCP data. It is 
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hypothesized that semi-diurnal and diurnal tides, winds and baroclinic tidal effects played 
a role in the forcing of these currents. 

B. PARO-SCIENTIFIC PRESSURE SENSOR 

A high precission Paro-Scientific pressure sensor was used to estimate surface 
elevation and to determine the mean depth of water. Surface elevation spectra were 
calculated by applying a linear wave transfer function to the pressure spectra. 



calculated from the pressure data averaged every five minutes. The data were acquired at 
2 Hz and with a thirty-minute ensemble. Fifty logarithmically-space frequency bins were 
calculated from 0.001 Hz to 0.2 Hz. Significant wave height was calculated from the 
surface elevation spectrum 



C. VIDEO CAMERA AND LASER LIGHT 

Underwater video cameras and structured light sources provide a means to 
observe sediment movement in near shore environments tmder conditions of good 
visibility. The video/laser system consists of two diode lasers and an imaging camera. 
The diode laser operates in the visible red at 685 nm wavelength with a power of 35 mW. 
The two lasers are housed in a rectangular pressure case 4.6 inches wide by 6 inches long. 
A cylindrical lens in front of each laser spreads the point sources into a 60° beam width, 
4 mm wide sheet of light. One laser sheet of light is directed in the cross-shore direction. 



=[ff,(/)]'s,(/), where /?,(/) = 



cosh(A:/z) 



. The depth of water was 



cosh(A:(/i + z) 




(3) 
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while the second is directed 45® to the cross-shore direction. The plane of light intersects 
the bed producing a bright stripe in the video image (see Figure 6a). 

A Cohu 4910, RS 170 monochrome camera is connected to one channel of the 
eight-node power and data controller at the MISO site, which provides serial data 
communication and a video bandwidth channel to shore. This underwater multi-node 
controller provides power to each node and converts the electrical video signal into a 
fiber optic signal for transmission 600 meters to the Marine Operations Building (MOB) 
on the shore. The electrical video signal is converted to a fiber optic signal to transmit the 
long distance while retaining low noise and high bandwidth. At MOB, the signal is 
converted back to a RS 170 signal and digitized by a Haupage fi:ume grabbing board in a 
Pentium III PC running the Linux operating system. The video firames are sampled at 10 
Hz. Each pixel in the images is processed to determine a mean, max, min, standard 
deviation and raw image stored over each five-minute ensemble. Initially the pixel size is 
640 by 480, but is reduced to 320 by 240 in order to conserve storage space and reduce 
noise. 

D. CALIBRATION 

The first step to extract meaningful position data from video images is to calibrate 
the field of view. To calibrate for underwater lens distortion, the camera was suspended 
in a laboratory water tank and exact measurements of distance from camera to a test grid 
were taken. To determine angle versus distance, calibration coefficients were computed 
in the cross-shore axis [cal(x)J and the ripple height axis [cal(z)] as shown in Figure 4. 
The calibrations from pixel space to angle for horizontal angle {fi) and vertical angle 
(«) is approximately given by 
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= [p/x(x) - 160] * cal{x) and a - [pixi^z) - 120] * cal{z) . 

The center pixel in each image was located at x=160 and z= 120 for the 320 by 240 pixel 
sized image. Following the geometry in Figure 5, distance from the camera to the bottom 
is determined by z = d* t 2 ca.{pff + a) , where d is the horizontal distance from the 
structured light to the focal point of the camera and off is the angle between the center 
pixel of the camera and the horizontal frame. For the data considered here, off is equal to 
47.48 degrees and d = 1 1 1 .76 cm. The cross-shore distance is estimated by 



x = - 






cos(<#) 

Since the instrument frame was tilted 12 degrees in the offshore direction during 
these observations, the raw image had to be rotated. Once rotated, the image was 
smoothed by filtering the data using a fourth order low-pass filter with a 2 cm cut off to 
produce a clean line representing relative ripple height in the cross-shore direction, with 
positive values of the cross-section signifying the offshore direction as displayed in 
Figure 6. The low-pass filter had negligible effect on measured slopes. 

During strong wave forcing, turbid conditions existed and video image resolution 
was significantly impaired. However, most of the time visibility was high, and the 
camera provided a high-resolution, two-dimensional time series of the cross-shore bed 
height over a 54 cm cross-shore distance with 1 .7 mm/pixel resolution. 
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III. RESULTS 



Wave induced ripple formation is related to the ability of waves to move 
sediment, which can be described by the Shields parameter. The Shields parameter is the 
ratio of the force mobilizing the sediments, measured by the bottom stress {t ^ ), to the 
restoring force of gravity 

where s is the ratio of sediment grain density to water density and is equal to 2.65 for 
quartz sand, g is the force of gravity, and is the mean grain size. Bottom stress can 
be divided into form drag and skin friction, where skin friction is primarily responsible 
for sediment motion and includes the effects of both waves and currents. Assuming a 
classic square law friction factor, and partitioning the stress into wave and current 
contributions, the equation for bed stress is: 

where is the rms wave velocity at the bed, U, V are the cross-shore and alongshore 
mean current velocities, and p is water density. The wave friction factor, 

= exp[5.213(— ^ — ^^0.194 _ 5 977 J jg based on Swart’s (1974) formula using the value 
A 

l.Sd^Q for the grain roughness of a flat bed (Fredsoe and Deigaard, 1992), where A is 
the orbital excursion radius. The current friction factor is /. = .003 (Amos et al., 1988). 
The initiation of sediment movement occurs when the grain roughness shields parameter 



11 



is greater than the critical shields parameter 6 ^^, = .04 (Grant and Madsen 1976), which 
is consistent with combined mean and oscillatory flows (Amos et al., 1988). 

The evolution of ripple patterns is displayed by examining sequences of the 5- 
minute mean images. Computer animated movies of the mean images provide a detailed 
view of these processes (http:/ /www.oc.nps .naw.mil/~stanton/ miso/dixon). The 22-day 
period of observation from yearday 245 through yearday 267 coincided with two wave 
events, two to three days in duration (Figure 7a). A wave event also preceded yearday 
245, which is included in the wave climatology. The wave-induced cross-shore bed 
velocity (Figure 7b) reflects the changes in wave height. The mean currents (Figure 7c) 
were weak with alongshore currents (<0.1m/s) stronger than the cross-shore currents 
(<0.05m/s). Diurnal variation in these currents (either due to tide or diurnal sea breeze) is 
evident. The combination of larger wave events succeeded by relatively small wave 
events created highly variable bedform evolution. 

Ripple wavelengths and migration are summarized in Figure 8a, which shows 
ripple height (grayscale intensity) versus cross-shore distance and time. The dotted lines 
are the crests of two ripples determined by a change in slope from positive to negative. 
Migration rates are determined from the slope of the crest line. The table describes ripple 
crest migration rates. 



Yearday 


246-248 


248.5-249.4 


254-257 


259.5-262 


264-265.3 


Crest# l(cm/day) 


6.84 


-5.6 


1.22 


6.3 


5.04 


Crest#2(cm/day) 


10.35 


-6.23/7.07 


1.89 
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Significant migration occurs between 246-248, 248.5-249.4 and 259.5-262. On yearday 
248.5, a bifurcation occurs as crest #2 separates, leaving one crest that continues to 
migrate offshore, while a second crest migrates onshore. On yearday 265, it is interesting 
to note that while forcing is low < 0 ^^, ), crest #1 migrates offshore fairly swiftly but 
crest #2 does not move. It is not understood what causes these differences in migration 
directions, however it is hypothesized that the crests are associated with a more three 
dimensional bedform. 

The Shields parameter was calculated using the measured BADCP velocities 
partitioned into wave mean square velocities (Equation 1) and mean current velocities. 
The wave velocities were approximately an order of magnitude greater than the mean 
current velocities. In addition, the wave fiiction factor is an order magnitude greater than 
the current fiiction factor. Hence, the contribution to the Shields parameter by the mean 
current was negligible (2 percent on average). 

Marginally active ripple patterns and evolution were seen throughout the period 
of observation despite the fact that the grain roughness Shields parameter (Equation 4) 
often dips below the critical Shields parameter = .04 (Figure 8 b). The video imaging 
system did not produce bedform images during daylight hours because the solar 
illumination usually swamped the structured laser light sheet. Therefore, each day is 
missing approximately 7 hours of images. Images are also missing during a large wave 
event between yeardays 249.5 and 252, because turbidity is strong during high sediment 
suspension. 

Ripple length A , is the distance fi:om crest to crest and ripple height, 7 , is the 

difference between the elevations of the crest and succeeding trough. A maximum angle 
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{<!>) of the ripple is calculated as the maximum angle occurring between a crest and 
succeeding trough, while the slope angle is the average angle between the crest and 

the difference 

in tan“' Ax = 0.8 cm being the image resolution. The slope angles were 

found not to be sensitive to the choice of Ax values, which was varied between 0.8 and 
2.4 cm. 

Images of the evolving ripple geometries are shown in Figure 9, where 3-4 day 
time series of data are plotted against relative ripple height and cross-shore distance. 
Observations begin with relic ripples during yearday 245, which are left over from the 
large wave days between yeardays 242-244. Relic ripples have smaller slope angles and 
more rounded crest than active vortex ripples. By yearday 246, 0 exceeds 0 ^^, and 
wave ripples begin to form. These are pointed, active vortex ripples, which shed vortices 
regularly. The Shields parameter remains above the 0 ^.^, threshold between yeardays 246 
and 249, during which time ripple migration is weakly offshore; however, due to the lack 
of bed profiles during these days, the ripple migration is difficult to resolve. The ripples 
are moderately asymmetric containing a steeper slope on the leeward (off-shore) side. 
During yearday 249, strong wave forcing caused strong sediment suspension, which 
blocked the vision of the video camera so that bed profiles were not recovered between 
yeardays 249.5 and 252. During this two and a half-day period, the Shields parameter 
reached a maximum over the twenty-two day time period. It cannot be determined if the 
ripples migrated onshore or offshore during this period of time 



succeeding trough,^ . Maximum angle (^) was determined by computing 
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On yearday 252, the waves calmed with^ below until yearday 254 leaving 
relic ripples, which significantly increase in size as only two crests dominate over the 57 
cm cross section. The ripple length of 30 cm between these successive crests is the 
longest during the entire period. Between yearday 252-254, the ripples are inactive 
representing relic ripples. From yearday 254-256, 6 is above the threshold value and 

bedforms begin to move slowly, migrating at a rate of 4-5 cm/day. This well organized 
long crested wave ripple is sustained as the bedform migrates seaward. 

Ripple migration ceases on yearday 256.5 as 0 is below and relic ripples 
prevail. The long crested wave ripple remains intact until yearday 258.5 when the second 
largest wave day begins with 9 well above 6 ^^, , and the ripples are active. Vortex ripples 
occur during the larger wave days, and their geometry is highly variable. Ripple 
migration begins at the same time that the large swell arrives (yearday 259). Migration 
rates are equal to that of the previous period at 5 cm/day. The long crested ripples that 
dominated earlier give way to shorter bifurcated three-dimensional ripples during 
yeardays 261 -263 (Figures 10a and 10b). Two ripple crests dominated the bed between 
252-259. The wave event (yearday 259) stirs the bed producing the bifurcation and three- 
dimensional ripples that are seen on yearday 261. Between yearday 261 and 265, wave 
energy is weak, and ripples are less active but continue to migrate a few cm/day. Bed 
movement slowly decreases and ripple migration comes to a halt on yearday 264. There 
is a small burst of wave energy on yearday 265, and ripples are highly active. From 
yeardays 266-267, wave energy is low and bedform variability is limited. The three- 
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dimensional bifurcated ripple gives way to symmetric gently sloping ripples on yearday 
267. 

In summary, the observation period is marked by three wave events. The first 
occurs prior to the start of ripple observation, the second wave event (249-252) is too 
strong to obtain bedform images, and bedform geometry is highly variable during the 
third wave event. It is during these larger wave events that ripple migration rates are the 
highest while ripples are strongly asymmetric with steeper slopes on the leeward side of 
the crest. Immediately following the larger wave events, relic ripples prevail as^ is 
below 0 ^, . Once wave energy increases slightly, ripples begin to migrate, although their 
geometry does not significantly change. The migrating vortex ripples yield to relic 
ripples as wave energy decreases and the pattern starts over. 
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IV. DISCUSSION 



Ripple steepness is limited by the angle of repose of the wet sand, where ^ is the 
local maximum angle between crest and succeeding trough. Assuming the ripple is a 
triangle, Nielsen (1981) suggests that the maximum steepness a bedform can obtain 

(^) = 0.5 tan ^ . If the trough had the shape of a parabola, the limit would be 

{rj! X)- 0.25 tan ^ as displayed in Figure 1 1 . Most of the bed slope data from the 
observations lies within the interval 0.25 tan^ < (r;/X) < 0.5tan^ (Figure 12). Nielsen 
(1979) fotmd that ripples maintain their maximum steepness when 0 < 0.2 but, when 
^>0.2, the vortex can no longer maintain the maximum steepness against the erosion of 
the crest. The result is that 7//X decreases as water particle velocity increases, and at 
stronger forcing ripples disappear completely. During the observation period, 0 peaked 
at 0.12 and we do not see ///A decreasing as a function of water particle velocity (Figure 
13a). Ripple steepness varied between 0.04 and 0.26 at low values of 9 <6 era > indicative 
of relic ripples. 

The ripple length has been suggested to be related to the orbital diameter (2A) of 
the wave-induced fluid motion near the bed, X = \.5A (Nielsen 1992). However, we find 

significant variation of^ over the range 9 < 0.05 , with ^ decreasing to an 

approximate value of 1.3 (Figure 13b). One has to be careful in the interpretation of the 

apparent decrease in ^ with increasing mobility number as 9 • 
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Therefore, both the ordinate and abscissa contain the same variable. Similar plots have 
been widely reported in the literature (Nielsen, 1992). For the larger value of 



approaches 1.3 independent of 6 in the range of ^>0.1 1 up to the maximum value of 



Clifton and Dingier (1984) summarized several data sets firom laboratory and field 
studies and suggested that that ripple wavelength scales linearly with orbital diameter for 

shorter period waves until orbital diameter/grain size ( ) reaches a value of 2000. 
These are described as orbital or vortex ripples (Bagnold, 1946). In the range of 



50 

go through a transitional stage, where ripples are much shorter and are referred to as 
suborbital ripples. For values of ^ , >5000, when orbital diameter is larger, ripple 



wavelength scales directly with grain size. The bedform environment was in a state of 
flux during the observation period. Orbital and suborbital ripples were the dominant 
ripple pattern (Figure 14). However, when wave forcing was small, relic ripples were 
present. The period starts with remnant or relic ripples, which are left over fi'om the 



fell between 2000 and 5000 on yearday 244. Ripples tend to orbital ripples on yearday 
246 as ripple length scales approach orbital excursion. The bed is in constant motion 
with ripple steepness peaking on yearday 246 and the steepness of these vortex ripples 
approach their angle of repose. When the largest waves of the observation period arrive 




0 . 21 . 




<5000, for longer period waves with increased orbital excursion, the ripples 




larger waves on yearday 243-244. These appear to be suborbital relic ripples as 2 
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on yearday 249.4, no images are retrieved, but it can be hypothesized that these are 
highly active suborbital ripples. When the waves subsided, suborbital relic ripples were 
left. The strong wave forcing changed the bedforms to two ripple crest spaced 30 cm 
apart. When the next large wave event occurs on yearday 259, ripple length scales with 
orbital excursion and orbital ripples are formed. Once the forcing ceases, orbital relic 
ripples remain and continue through to the end of the observation period. 

Ripples become asymmetric when suspended sediments are lifted and carried 
from the upstream side trough and settle on downstream side crest. In the following 
discussion, upstream is the side of the ripple where the sediment is moving from and 
downstream is the location where the sediment is moving. The building up of the 
downstream side crest creates a steeper slope on the downstream side and a shallower 
slope on the upstream side. As this process continues, the downstream side becomes 
steeper and the upstream side becomes shallower. Active asymmetric ripples migrate in 
the direction of the steeper slope, while active symmetric ripples may not change their 
position. For example, during yearday 254, the leeward slope is twice as steep as the 
windward slope indicative that the ripples are migrating in the offshore direction (see 
histograms of slope angle for selected ripple profiles in Figure 15). 

Mean slope angles were calculated by subtracting the mean slopes of the offshore 
and onshore sides of the ripple to determine ripple symmetry. Large positive values 
represent asymmetric ripples with a steeper slope on the offshore side, indicative of 
ripples migrating offshore (negative slopes indicate onshore migration) while small 
values correspond to symmetric ripples (Figure 16). Asymmetric ripples were formed 
following the wave events and slowly eroded with weaker forcing. The forcing by larger 
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waves on yeardays 243-244 appears to have resulted in the observed asymmetric ripples 
on yeardays 245-246. These bedforms became symmetric as the wave forcing subsided. 
The larger waves starting on yearday 248 had large onshore, negative skewness (see 
Figure 17) coinciding with onshore migration (see Figure 8a) and ripple asymmetry 
indicative of onshore transport. From yeardays 250-253, 0 decreased and skewness 
changed from onshore to offshore. When the images resumed, ripples migrated onshore 
(yeardays 253-257) and were asymmetric. As wave forcing decreased, migration rates 
slowed and ripples became symmetric. During the last wave event (yeardays 259-262) 
migration rates were strong. Ripples became asymmetric for a short time on yearday 263 
and tended toward symmetry as the period of observation concluded. 

Mean and skewness velocities were calculated to assess if they could explain 
offshore ripple migration. Skewness of low-passed and high-passed filtered velocities 
(0.05 Hz cut-off), most of the time showed no preferred direction (Figure 17). The 
exceptional times when onshore, negative skewness exceeded 0.1 occurred on yeardays 
250 and 265 when migration was onshore. Cross-shore mean current were weak and 
fluctuated between the offshore and onshore direction, again suggesting no preferred 
direction (see Figure 7c). 

Bound long waves give an offshore contribution to sediment transport and can 
result in offshore ripple migration as schematically shown in Figure 18. The presence of 
bound long waves was measured by correlating the infragravity wave band (obtained by 
low passing the cross-shore velocity below 0.05 Hz) with the envelope of the sea-swell 
band of frequencies (Hanes 1991) as shown in Figure 19a. The envelope of the sea-swell 
band is obtained by low pass filtering (below 0.05 Hz) the time series. Positive 



20 



correlation suggests bound long waves, (coordinate system is u positive offshore). The 
correlations are significantly different than zero at the 95% confidence level (assuming 
Guassian velocities) for-0.0167< r <0.0167. The r values are persistently positive and 
mostly significant indicating bound long waves were present 

A measure of sediment transport advected by bound long waves is given by 

u,p * where sediment transport is parameterized as proportional to (w^)/^ (Figure 

19b). Transport was mostly offshore in agreement with ripple migration direction. 

Four offshore ripple migration events occurred between 246-248, 254-257, 259.5- 
262, and 264-265.3 as shown in Figure 9. The first event yielded the swiftest migration 
rates. Orbital forcing was strong enough for^^^;, threshold to be surpassed. Two ripples 
dominated at the start of yearday 246 and the ripple #2 migrated at 10.35 cm/day while 
ripple #1 migrated at 6.84 cm/day. During this time frame, waves were groupy and bound 
long waves were present. The migration rate was low during the second event as 6^^, was 
exceeded only during yearday 255 by 0.02. These ripples migrated offshore at 1 .5 
cm/day. Bound long waves were present and ripples slowly migrated offshore. Yeardays 
259.5 and 262 provide good information because although the forcing was strong, video 
images were obtained. The presence of bound long waves caused the ripples to migrate 
offshore at 6.3 cm/day. 

There is a question of what causes the final event to migrate offshore at 5.0 
cm/sec when 6 < . This brings into question the validity of 0^^, =0.04 or how this 

value is applied. By viewing the video showing ripple forms migrating along with 
calculated Shields parameter, it is qualitatively observed that the bed only migrated when 
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6 > with exception of yearday 264. The basis for ^^,,,=0.04 are laboratory 
experiments using monochromatic waves with constant amplitude (Grant and Madsen, 
1976). The wave velocity is a bulk rms value for a 30 minute ensemble and does not 
accurately capture the large range of forcing between the peak and minima within a wave 
group. For very narrow banded, groupy waves in nature, the maximum wave velocity 
could be more two-five times the monochromatic value for the same Shields Parameter. 
Therefore, it is hypothesized 6 ^^, may be less than 0.4 for waves in nature. 
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V. CONCLUSION 



The analysis of bedfonn evolution using an underwater video camera in 
conjunction with a laser light has been demonstrated. High-resolution images were 
obtained and analyzed in relation to wave and current forcing. Ripple geometry and 
migration rates were obtained. The follovving conclusions are made. 

Ripple steepness varied for values of 6 , and no relationship was found. Relative 

ripple length (^ ) decreases as a function of the Shields parameter and was 

approximately constant at 1 .3 in the range of 0.07 <0 < 0.21 . Ripples were steeper on 
the downstream or offshore side. Asymmetric ripples were left after large wave events 
and the tended toward symmetry with weaker forcing. Orbital ripples were the dominant 
bedform found when the wave and current forcing was strong. However, when the waves 
relaxed relic ripples prevailed. Forcing was never strong enough such that orbital 
excursion could reach a value to create anorbital ripples, however, suborbital ripples were 
formed during the strongest forcing event. 

Migration of ripples may be the dominant mode of cross-shore sediment transport 
at the MISO site. Measured migration rates occurred in four events (246-248, 254-257, 
259.5-262, and 264-265.3), and all ripples actively migrated offshore. During these 
events, B exceeded the critical shields parameter on all days except during yearday 264- 

265.3. Skewness of orbital velocity (both high and low passed) mostly showed no 
preferred direction to explain offshore ripple migration. Significant positive correlation 
was obtained between wave groupiness of the narrow banded swell, (measured as the 



23 



envelope of high passed velocity signal), and bound long waves, which may account for 
the offshore migration during the first three events. 

The use of video imaging to accurately map the bedform evolution was 
successful, however when large waves caused strong sediment mixing, no video images 
were obtained. Hence ripple information was limited to low wave events. 
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Figure 1. Instrument frame (viewed from sea to shore) at MISO array. ADCP was used to 
measure current velocities. The PARO was used to calculate wave height and water 
depth. The camera and structured light provided high resolution video images of 
bedforms. 
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Figure 2a. Location of MISO instrument array. The star represents the location of the 
MISO array. Figure 2b. Regional map of central coast of California showing that the 
Point Pinos headland protects the MISO location and causes waves to be narrow-banded 
in frequency and direction owing to strong refraction. 
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Sand Fraction 



Sand Sediment Distribution at MISO Site - 




Figure 3. Sand Grain Size Distribution with median grain size )=.125 mm 
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Actual Angle (radians) Actual Theta (radians) 



Actual Angle (Radians) vs Horizonteil Distance (Pixels) 





Vertical Distance (Pixels) 



Figure 4a. Horizontal calibration Figure 4b. Vertical calibration. The slope of the line is 
the calibration coefficient. 
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d=l 11.76 





Y=d * tan(B )/cos(off) 
B=(pix(x)-160)*cal(x) 




Figure 5. Camera calibration showing horizontal and vertical cross-sections. Alpha and 
B are angles in the vertical (z) and cross-shore (y) directions. 
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Figure 6a. Raw image taken from video camera. Figure 6b. Conversion of video image 
in pixel space to raw image in real space (cm). Figure 6c. Image was rotated because the 
instrument frame was tilted 12 degrees offshore. Figure 6d, The image was then 
smoothed by applying a low pass filter with a 3cm scale. The filter did not affect the 
image resolution which was 1.7mm/pixel. 
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Figure 7a. Significant Wave Height ) for yearday 243-267 of 30-minute averages. 
Figure 7b. Rms orbital velocity at the bed ( ) of 30-minute averages. Figure 7c. 
Mean currents (cross-shore and alongshore of 1-hour averages). 
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Ripple Migration 
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Figure 8a. Time series of Relative ripple height displayed in grayscale (white is large, 
dark is small) vs. cross-shore distance (positive values are offshore) and time. Ripple 
Crest #1 is the lower ripple crest and Ripple Crest #2 is the upper ripple crest. Greatest 
migration rates occur between 244-248. Figure 8b. Grain roughness Shields parameter 
calculated from a rms orbital velocity averaged every 30 minutes. 
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Figure 9. Time evolution of ripple geometries displayed in 3-4 day periods of relative 
ripple height (grayscale) vs. cross-shore distance (positive is offshore). 
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Figure 10a. Raw image of long crested ripples on yearday 253 (9/10/99) evolve into 
shorter bifurcated ripples (Figure 10b) on yearday 261 (9/18/99). The long crested two 
are left after the strong forcing on yearday 249-252. These give way after the high waves 
on yearday 259-261 to a short crested bifurcation. 
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Figure 1 1. With maximum slope on the ripple profile is equal to tan ^ , the height to 
length ratio must for geometric reasons be 0.5 tan (p for a triangle ripple and 0.25 tan <l> 
for a parabolic ripple (after Nielsen 1992). 
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Ripple Geometry 




Figure 12. Time series of ripple steepness/tan 0 . Bed slope geometry data generally falls 



between 0.25< 




tan<z> 



/X 

tan^ 



<0.25 when is not exceeded. 
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Figure 13a. Ripple steepness is highly variable. Figure 13b. Relative ripple length 
parameterized by orbital excursion (plotted only when 9 > 9 ^^, ) decreases as forcing and 
orbital excursion increases. The dashed line corresponds to A = 1.3 A . 
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Figure 14a. Ripple Geometry (Clifton and Dingier, 1976) classified into orbital. 



(2^ ) < 2000, suborbital, 2000 < 2^ < 5000, and, anorbital, 2^ 



>5000. Figure 
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14b. Time series of A and A . Only values of A have been plotted when 9 > ^cn,- 
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Figure 15. Left panel show evolution of selected images of cross-shore distance (positive 
is offshore) vs relative ripple height. Right panel shows histogram of slope angle with 
mean slope angle displayed above histogram. Mean slope angles were calculated by 
averaging the slope angle between crest and succeeding trough. The negative slopes 
angle are the mean slopes on the off shore side, while the positive slope angle are the 
slope angles on the onshore side. 
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Figure 16a. Ripple Asymmetry: Large positive values show steeper sloping ripples on the 
off-shore side indicative of offshore migration while small values show symmetric 
ripples. Negative values show steeper sloping ripples on the on-shore side. Asymmetric 
ripples were left after strong wave forcing (high value of Shields Parameter, Figure 16b) 
and tended toward symmetry as the forcing subsided. 
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Skewness of Surface Gravity Waves 




yearday 



Figure 17. Skewness of cross-shore orbital velocity from surface gravity waves. A high 
pass filter at 0.05 Hz was applied to the cross-shore velocity time series and skewness 
values were calculated over one-hour samples. 
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BOUND LONG WAVE TEffimY 





Figure 18 . Change in momentum (high waves) is balanced by set down of forced long 
wave and small waves is balanced by a relative set-up. The forced long wave gives 
reverse flow under high wave and forward flow under small waves. Because more 
sediment is suspended during the larger waves. Net transport is revered or offshore. 
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Figure 19a. Correlation Coefficient determined by correlating the infragravity wave band 
(obtained by low pass filtering the cross-shore velocity at 0.05 Hz) with the envelope of 
the sea-swell band of frequencies. The dashed lines represent the significant correlation. 
Figure 19b. Correlation as calculated in (a), but is not normalized. Correlation is 
determined by u,p * {u^)ip ■ 



43 



TfflS PAGE INTENTIONALLY LEFT BLANK 



44 



LIST OF REFERENCES 



Amos. C.L., Bowen, A.J., Huntley, D.A., and Lewis, C.F.M., Ripple Generation Under 
the Combined Influences of Waves and Currents on the Canadian Continental Shelf. 
Cont. Shelf Res., 8, 1129-1153, 1988. 

Andersen, K.H., The Dynamics of Ripples Beneath Surface Waves and Topics in Shell 
Models of Turbulence, PH. D. dissertation. University of Copenhagen, 142, pp., 1999. 

Bagnold, R.A., Motions of Waves in Shallow Water; Interaction Between Waves and 
Sand Bottoms, Proc. R. Soc. London, Ser. A, 187, 1-15, 1946. 

Boyd, R., D.L. Forbes, and D.E. Heffler, Time-Sequence Observations of Wave-Formed 
Sand Ripples on an Ocean Shoreface, Sediment ology, 35, 449-464, 1988. 

Clifton, H.E., and J.R. Dingier, Wave-formed Sedimentary Structures and 
Paleonenvironmental Reconstruction, Mar. Geol., 60, 165-198, 1984. 

Crawford, A.M., and A.E. Hay, A Simple System for Laser-illuminated Video Imaging 
of Sediment Suspension and Bed Topography, lEEEJ. of Ocean. Eng., 23(1), 12-19, 

1998. 

Deigaard, R., J.B. Jakobsen, and J. Fredsoe, Net Sediment Transport Under Wave Groups 
and Bound Long Waves, J. ofGeophys. Res., 104 (C6), 13559-13575, 1999. 

Engelimd, F., and J. Fredsoe, Sediment Ripples and Dimes, Annu. Rev. Fluid Mech., 14, 
13-37, 1982. 

Fredsoe, J., and R. Deigaard, Mechanics of Coastal Sediment Transport, World Sci., 
River Edge, N.J., 1992. 

Grant, W.D., and O.S. Madsen, Quantitative Description of Sediment Transport by 
Waves, Proc. 15'^ Int. Conf On Coastal Engineering, Honolulu ASCE, 1093-1112, 1976. 

Hanes, D.M., Suspension of Sand Due to Wave Groups, J. Geophys, Res., 96 (C5), 8911- 
8915, 1991. 

Hay, A.E., and D.J. Wilson, Rotary Sidescan Images of Nearshore Bedform Evolution 
During a Storm, Mar. Geol., 119, 57-65, 1994. 

Inman, D.L., and A.J. Bowen, Flume Experiments on Sand Transport by Waves and 
Currents, in Coastal Engineering, Proceedings of the 8‘^ Conference, pp. 137-150, Am. 
Soc. Civ. Eng., New York, 1963. 



45 



Longuet-Higgins, M.S., and R. W. Stewart, Radiation Stress in Water Waves: A Physical 
Discussion with Applications, Deep Sea Res., Ill, 529-562, 1 964. 



Nielsen, P., Some Basic Concepts of Wave Sediment Transport, Ser. Paper 20 Inst 
Hydrodyn Hydraul Eng, Tech Univ Denmark, 160 pp., 1979. 

Nielsen, P., Dynamics and Geometry of Wave Generated Ripples, J. Geophys. Res., 86, 
C7, 6467-6472, 1981. 

Nielsen, P. Coastal bottom Boundary Layers and Sediment Transport, World Scientific, 
1992. 

Shi, N. C., and L.H. Larsen, Reverse Sediment Transport Induced by Amplitude 
Modulated Waves, Mar Geol., 54, 181-200, 1984. 

Swart, D. H., Offshore Sediment Transport and Equilibrium Beach Profiles, Delft Hydr 
LabPubl,U\,\91A. 

Stanton, T.P., WWW.oc.nps.navy.mil/~stanton/miso 

Tanaka. H., and Shuto, N., Sand Movement Due to Wave-current Combined Motion. 
Coast. Engrg. in Japan, Tokyo, 27, 179-191, 1984. 

Traykovski, Peter, A.E. Hay, J.D. Irish, and J.F. Lynch, Geometry, Migration and 
Evolution of Wave Orbital Ripples at LEO-15, J. Geophys. Res., 104, Cl, 1505-1524, 
1999. 



46 



INITIAL DISTRIBUTION LIST 



No. Copies 

1 . Defense Technical Information Center 2 

8725 John J. Kingman Rd, STE 0944 

Ft Belvoir, VA 22060-6218 

2. Dudley Knox Library 2 

Naval Postgraduate School 

411 Dyer Rd 

Monterey, CA 93943-5101 

3. Chairman (Code OC/Gd) 1 

Department of Oceanography 

Naval Postgraduate School 
Monterey, CA 93943-5122 

4. Chairman (Code MR/Wx) 1 

Department of Meteorology 

Naval Postgraduate School 
Monterey, CA 93943-5114 

5. Dr. Edward B. Thornton, (Code OC/Tm) 2 

Department of Oceanography 

Naval Postgraduate School 
Monterey, CA 93943-5122 

6. Professor Tim Stanton, (Code OC/St) 2 

Department of Oceanography 

Naval Postgraduate School 
Monterey, CA 93943-5122 

7. Damon Dixon , 4 

704 San Gabriel Place 

Chula Vista, CA 91 914 

8. Lori Boettler 1 

16 Greenwood Court 

Orinda, CA 94563 

9. Cameron Hansen 1 

7693 Summerhill Court 

Salt Lake City, UT 84121 

47 



72 290NP6 



6/02 



TH 

22S27-: 



’3013 

■200 nlb 



